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Abstract 

The polarized electromagnetic wave in cylindrical fibers is 
described by a system of four component vector nonlinear 
Schrodinger equations. We have investigated analytically 
and numerically the modulational instability of this system. 
Several types of modulational instability are shown to exist 
depending on the sign of the dispersion coefficients. The 
results indicate that the group- velocity mismatch and the 
four wave mixing terms influence considerably the shape, 
position and especially the number of the spectral regions of 
the gain spectra of the modulation instability. Numerical 
simulations of the four component vector nonlinear 
Schrodinger equation show that the phenomenon of 
modulational instability gives rise to periodic pulse arrays of 
waves train. The presence of the four wave mixing terms 
influences considerably the evolution of the periodic pulse. 
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Introduction 

Vector solitons in multiple waveguides have drawn 
tremendous attention recently. While solitons can 
propagate over a long distance from a balance of 
dispersion and nonlinearity, coupling in two or more 
waveguides [1] permits additional flexibility which 
generates a rich variety of phenomena. The study of 
the dynamics of nonlinear waves is central to 
understand the behaviour of a wide range of physical 
systems, taking reference to [2], for example. In recent 
years, an important theme of much of this work has 



been the study of modulational instabilities (MI) of 
plane waves. MI occurs as a result of the interaction 
between the linear dispersive effect and the nonlinear 
effect, which is characterized by instability 
experienced by a continuous wave (CW) when it 
propagates inside a dispersive nonlinear system 
together with a weak noise. The perturbed wave 
experiences experimental growth and this leads to the 
breakup of the temporal profile of a wave into a train 
of ultrashort pulses [3, 4]. Since this disintegration 
typically occurs in the same parameter region where 
bright solitons are observed, MI is considered to some 
extent, a precursor to soliton formation in nonlinear 
systems [3]. The MI gain has been extensively studied 
for the scalar nonlinear Schrodinger type equations 
(NLSEs) [4]. The experimental evidence was then 
given that MI can be observed in a non-birefringent 
fiber with anomalous group velocity dispersion (GVD) 
[5]. This result has been confirmed subsequently by 
other observations of MI in the anomalous dispersion 
regime [6]-[7]. Vector MI can occur in isotropic fibers 
such that the gain spectrum depends on the 
polarization state of the input CW beam. It has been 
known, however, since the early studies by Bespalov 
and Talanov [8], and Berkhoer and Zakharov [9], that 
an incoherent (i.e., purely intensity dependent) 
couplingbetween two NLSE leads to extension of the 
instability domain in the normal dispersion regime. 
Unfortunately, it is very difficult to manufacture 
birefringence-free fibers. As an alternative, MI was 
observed in a bimodal fiber in which the input beam 



1 



www.seipub.or^rap 



Review of Applied Physics (RAP) Volume 2 Issue 1, March 2013 



excited two fiber modes with nearly equal power 
levels, and the two modes had the same group 
velocity [10]. However in a multi-mode fiber the 
velocity difference (walk-off) has to be taken into 
account, because different modes have substantially 
different group velocities. Previous studies have 
shown that the effect of walk-off among optical pulses 
can considerably modify the gain spectra of cross- 
phase modulation (XPM) induced MI [11]. So far, a 
different situation that is of considerable interest in 
nonlinear optics is the coherent coupling between two 
polarizations in a weakly anisotropic media. Coherent 
coupling can arise in the polarization of light beams, 
as the relative phase factors of the interacting electric 
fields there play a crucial role in the form of four-wave 
mixing (FWM) process. FWM in optical fibers can be 
both harmful and beneficial depending on the type of 
application. It can induce crosstalk in wavelength 
division multiplexing communication systems and 
limit the performance of such systems. The FWM- 
induced crosstalk can be minimized in practice 
through dispersion management, a technique in which 
the dispersion of each fiber section is made large 
enough that the FWM process is not phase matched 
through-out the link length [12]. At the same time, 
FWM is useful for variety of applications [13]. FWM 
can be useful for applications such as optical sampling, 
channel demultiplexing, pulse generation, and high- 
speed optical switching [14]. It can also be used for 
reducing quantum noise through squeezing, and for 
generating photon pairs that are quantum-correlated. 
The resulting model is a system of coherently coupled 
NLSE (CNLSE). In the light of the above mentioned 
facts, it is worthwhile to study the phenomenon of MI 
in a fiber system with the FWMnonlinearity. Recently, 
MI in a system of coherently CNLSE has been 
investigated and it was found that the nonlinearity 
arising from phase coherence enlarges the range of 
XPM- induced MI [15]. Modulation instabilities arising 
from systems of four CNLSE have also been studied 
earlier in the literature [16]. It was found here that MI 
can appear even in the regime where the XPM 
coefficient is less than that of the self-phase 
modulation (SPM). The group-velocity mismatch 
(GVM) and the mixed interaction between different 
modes (FWM) have been neglected. Therefore, 
generally speaking, both of the walk-off effects and 
FWM terms should be taken into consideration in the 
study of XPM induced-MI in the four CNLSE. Our 
motivation in the present work is to study the MI of 
the four CNLSE, derived by Leble and Reichel [17]. 



The remainder of the article is arranged as follows. In 
section 2, a system of four CNLSE is presented. In 
section 3, we study the linear stability analysis of this 
model equation. The results of MI gain are 
investigated in section 4. Typical outcomes of the 
nonlinear development of the MI are reported in 
section 5 and section 6 concludes the paper. 

The model 

Wave propagation for pulses of polarized 
electromagnetic waves in cylindrical fibers is governed 
by a CNLSE [17] 

dX& + dX& . d 2 x£ r I + I 2 III 2 

lf =Q; o ± i^r + <7^-4 iKil +o 2 \x \\ + 
3 \X± | 2 + 4 \X+ | 2 + OsX^X* + o 6 x*xR]x± v 

dX& . dxA + d 2 x£ r I . .2 iii 2 

it =a nir +tfu-#-*N*?il +( ? 2 fel + 

2-2 — 
1^01 I + 1-^01 I + $5^01 -^01 + @6^0* ^01*]-^ir 0- ) 

Here, and are the pulses envelope in each 
polarization mode (two modes and each mode has two 
polarizations).^ and /?.* (j=0, 1) denote the GVM of 
pulses and the GVD, respectively. 1 and Q 1 are SPM 
coefficients. 2 , 3 , 4 , Q 2 , Q 3 , and Q 4 are the XPM 
coefficients. The nonlinear coefficient s0 5 , 6 ,Q 5 and Q 6 
describe mixed interaction between different modes 
and different polarization (FWM). In the limit of high 
birefringence, the FWM mixing terms in Eq. (1) can be 
neglected as they can be averaged to zero. In contrast, 
these coherent phase nonlinearities will be retained in 
the present work. Although solitons have been studied 
intensively [18-21], the focus here will be on the MI of 
plane waves of Eq. (1). Many existing works have 
assumed the low birefringence limit. The first objective 
here is to demonstrate that restoring this physics 
changes considerably the MI sidebands. Secondly, 
inclusion of phase dependent nonlinearity will be 
shown to extend the XPM- induced MI of the CNLSE. 

The linear stability analysis 

In order to investigate the evolution of weak 
perturbations along the propagation distance, we 
carry out a linear-stability analysis. Firstly, we look for 
CW solution of the system described by Eq. (1 ) given 
by: 

X*=>flj*exp®-t0,z), j=0,l (2) 

where P.* is the power of the plane wave and the real 
phase0 ; - are given by: 
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O = O 1 |P±| 2 +O 2 |P+| 2 + 3 \P*\ 2 + 4 |P+| 2 + Orjp + P* + 



Qi=Qim + Q2\Pii\ +Q 3 K \ + Q4\p + i \ + Q5JOT1 + 



QejPoiPoi. 



(3) 



In order to examine the linear-stability of the CW 
states, we introduce the perturbed fields of the form 



^•T=(fe ±+ a/l)exp(-i0 7 -z), j=0, 1 



(4) 



Where « ^jPj^ . By substituting Eq. (1) and 

retaining linear terms of ay^, we obtain four linearized 
equations: 

= a± "-f + l P ± ^ - i[ 0l P± (a± + a£) + 



°2j«t ( a + i + O + 3 J^iT^ ( a++a+*) + 
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(5) 



where * denote complex conjugate. We assume 
general solutions of the form 



<^i = w*exp(iOt) +K*exp(— iQt), j=0, 1 



(6) 



Inserting Eq. (5), one obtains a set of eight linear 
ordinary differential equations 



d[Y]_ 

dz 



i[M][Y] 



(7) 



where [M] is an 8x 8 matrix with 



[Y]=[uJ 1 ,Voi' u oi> v oi> u ti> v u> u ii> v \iY ■ The elements 
of [M] are given in appendix. The above 8x8 matrix is 
referred to as stability matrix which is used to study 
the stability of the system under consideration. MI 
occurs only when at least one of the eigenvalues of the 
stability 8x 8matrix possesses a nonzero and negative 



imaginary part, which results in an exponential 
growth of the amplitude with the perturbation. By 
requiring the determinant of the associated matrix in 
Eq. (7) to vanish, we can find the solvability condition 
for these equations, which amounts to a polynomial of 
degree eight. The MI phenomenon is measured by a 
power gain given by G = 2\lm(K)\ where K is the 
eigenvalue of MI that possesses the largest imaginary 
part [22]. 

The Modulational Instability Gain 

The MI gain will be calculated and discussed for 
different dispersion regions in the following. To 
simplify the discussion, we set the incident optical 
powers of the four optical waves as = P. We set 
also the GVM as = — tt£i= 5. We have also chosen 
the following physical parameters [23]: (3j^ = 
0.06ps 2 /m()=0, 1), 1 = Qj= 1, 2 = Q 2 = 2, 3 = 4 = 



Q 3= Q 4 = 2/3, 5 



Or, 



1/3, P = 1.5 W. The 



parameter 5 = 2—5 ps/m for a wavelength difference 5 
- lOnm. Several qualitatively different situations 
emerge depending on the signs of dispersion 
coefficients. 

Case 1: fa < 0,fa < O,/^ < 0,/?£i < 

In this case, all the four beams propagate in the 
anomalous GVD of the fiber. This case is of particular 
importance since each beam can develop MI even in 
the absence of the XPM coupling. The gain spectra of 
MI for different parameters 5 have been shown in FIG. 
1. We have obtained in this case in FIG. 1 (al) only one 
sideband. 




FIG. 1(A1) THE MI GAIN SPECTRA WITHOUT THE FWM 
NONLINEARITIES. THE FOUR BEAMS EXPERIENCES 
ANOMALOUS GVD. THE FIBER PARAMETERS ARE /? ± = 
-0.06Ps 2 /M /?* = - 0.06Ps 2 /M, t =&= 1, 2 = Q 2 = 2,0 3 = A = 
Q 3 = Q t = 2/3, 5 = 6 = Q 5 =Q 6 = 1/3, P= 1.5 W 

As 5 increases, the number of the spectral peak 
increases also as seen in FIG. 1 (a2), 
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FIG. I(a2) THEMI GAIN SPECTRA WITHOUT THE FWM 
NONLINEARITIES . THE FOUR BEAMS EXPERIENCES 
ANOMALOUS GVD. THE FD3ERPARAMETER5 ARE /?± = 
-0.06ps 2 /m, = - 0.06ps 2 /m, O t = 1, 2 = Q 2 = 2, 3 = 4 = 
Q 3 = (? 4 = 2/3,0 5 = 6 = Q 5 =Q 6 =1/3,P = 15W 

and on further increasing 6, we observe in FIG. 1 (a3) 
that the additional sidebands are suppressed. 




FIG. I(a3) THEMI GAIN SPECTRA WITHOUT THE FWM 
NONLINEARITIES. THE FOUR BEAMS EXPERIENCES 
ANOMALOUS GVD. THE FIBER PARAMETERS ARE /?* = 
-0.06ps 2 /m, /?± = - 0.06ps 2 /m, O l = 1, 2 = Q 2 = 2, 3 = 4 = 
Q 3 = Q 4 = 2/3, 5 = 6 = Q 5 =Q 6 =1/3,P = 15W 

For the large values of 6, we observe in FIG. 1 (a4) only 
one sideband. 




By comparison in FIG. 1 (b) where we have taken into 
account the FWM nonlinear coefficients, the gain 
spectrum takes on new characteristics and the number 
of the spectral peak are quite different from those of 
FIG. 1 (a). 




FIG. l(bl) THE MI GAIN SPECTRA IN THE CASE WHERE THE FWM 
NONLINEARITIES ARE TAKEN INTO ACCOUNT. THE FIBER 
PARAMETERS ARE /?± = -0.06ps 2 /nv/?u= - 0.06ps 2 /m, O l = 1, 
2 = Q 2 =2,0 3 =0 A =Q 3 = Q A =2/3,0 5 = 6 =Q S =Q 6 =1/3,P = 
1.5 W 




FIG. I(b2) THE MI GAIN SPECTRA IN THE CASE WHERE THE FWM 
NONLINEARITIES ARE TAKEN INTO ACCOUNT. THE FIBER 
PARAMETERS ARE = -O.OepsVm,/?^ - 0.06ps 2 /m, O l = 1, 



2 = Q 2 =2,0 3 =0 4 



= Q 3 = Q 4 =2/3,0 5 = 6 
1.5 W 



<2 6 =1/3,P = 




FIG. I(a4) THEMI GAIN SPECTRA WITHOUT THE FWM 
NONLrNEARITIES . THE FOUR BEAMS EXPERIENCES 
ANOMALOUS GVD. THE FIBER PARAMETERS ARE /?± = 
-0.06ps 2 /m, = - 0.06ps 2 /m, O t = Q t = 1, 2 = Q 2 = 2, 3 = 4 = 
Q 3 = Q 4 = 2/3,0 5 = 6 = Q S =Q 6 =1/3,P = 1.5W 



FIG. I(b3) THE MI GAIN SPECTRA IN THE CASE WHERE THE FWM 
NONLrNEARITIES ARE TAKEN INTO ACCOUNT. THE FIBER 
PARAMETERS ARE p± = -0.06ps 2 /m, Prr ~ 0.06ps 2 /m, O l = Q t = 1, 
2 = Q 2 =2,0 3 =0 A =Q 3 = Q A =2/3,0 5 = 6 = Q s = Q 6 = 1/3, P = 
1.5 W 
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FIG.l(b4)THE MI GAIN SPECTRA IN THE CASE WHERE THE FWM 
NONLINEARTHES ARE TAKEN INTO ACCOUNT. THE FIBER 
PARAMETERS ARE /?* = -0.06p> 2 /m, - 0.06ps 2 /m, o l = &= 1, 
2 = Q 2 =2,0, =0 4 = Q 3 = Q 4 = 2/3,0 5 = 6 =Q 5 = Q 6 =1/3,P = 

1.5 W 

Case 2: fa < 0,0^ < 0,^ > 0,/?^ > 

In this case, the first mode experiences anomalous 
GVD, while the second mode propagates in the 
normal GVD region. When 5 = 2, we observe in FIG. 2 
(al) only one spectral peak, 




FIG. 2(al) THEMI GAIN SPECTRA WITHOUT THE FWM 
NONLINEARTHES. THE FIRSTMODE EXPERIENCES 
ANOMALOUS GVD, WHILE THE SECOND MODE PROPAGATES 
IN THE NORMAL DBPERSION REGIME. THE FIBER 
PARAMETERS ARE /?* = -0.06ps 2 /m, = 0.06ps 2 /m. OTHER 
PARAMETERS ARE IDENTICAL TO THOSE USED FOR FIG. 1 




FIG. 2(a2) THEMI GAIN SPECTRA WITHOUT THE FWM 
NONLINEARTHES. THE FIRSTMODE EXPERIENCES 
ANOMALOUS GVD, WHILE THE SECOND MODE PROPAGATES 
IN THE NORMAL DBPERSION REGIME. THE FIBER 
PARAMETERS ARE = -0.06ps 2 /m, = 0.06ps 2 /m. OTHER 
PARAMETERS ARE IDENTICAL TO THOSE USED FOR FIG. 1 



but with the increase of the parameter 5, we observe in 
FIG. 2 (a2) that the position of the spectral peak varies, 
and in FIG. 2 (a3) that the number of the spectral 
region increases. 




FIG. 2(a3) THE MI GAIN SPECTRA WITHOUT THE FWM 
NONLrNEARITIES . THE FIRSTMODE EXPERIENCES 
ANOMALOUS GVD, WHILE THE SECOND MODE PROPAGATES 
TN THE NORMAL DBPERSION REGIME. THE FIBER 
PARAMETERS ARE = -0.06ps 2 /m, 0.06ps 2 /m. OTHER 
PARAMETERS ARE IDENTICAL TO THOSE USED FOR FIG. 1 

For 5> 5, we have seen that the spectral peaks are 
totally suppressed. When we take into account the 
FWM nonlinear coefficients, we observe in FIG. 2 (bl)- 
FIG. 2(b3), that these coefficients alter the number and 
the position of the spectral peaks. 




FIG. 2(bl) THE MIGAIN SPECTRA IN THE PRESENCE OF THE FWM 
NONLrNEARITIES. THE FIRSTMODE EXPERIENCES 
ANOMALOUS GVD, WHILE THE SECOND MODE PROPAGATES 
PJ THE NORMAL DBPERSION REGIME. THE FIBER 
PARAMETERS ARE = -0.06ps 2 /m, /?*= 0.06ps 2 /m. OTHER 
PARAMETERS ARE IDENTICAL TO THOSE USED FOR FIG. 1 

Case 3: > 0,/? "i > 0,/?^ < 0,/Jfj < 

In this case, the first mode experiences normal GVD, 
while the second mode propagates in the anomalous 
GVD region. The gain spectra of MI for different 
parameters 5 have been shown in FIG. (3). When 5=2, 
the gain spectrum in FIG. 3 (al) contains four spectral 
regions. 
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FIG. 2(b2) THE MIGAIN SPECTRA IN THE PRESENCE OF THE FWM 
NONLINEARTHES. THE FIRSTMODE EXPERIENCES 
ANOMALOUS GVD, WHILE THE SECOND MODE PROPAGATES 
IN THE NORMAL DBPERSION REGIME. THE FIBER 
PARAMETERS ARE = -0.06p> 2 /m, = 0.06ps 2 /m. OTHER 
PARAMETERS ARE IDENTICAL TO THOSE USED FOR FIG. 1 




FIG. 2(b3) THE MIG AIN SPECTRA IN THE PRESENCE OF THE FWM 
NONLINEAR1TIES. THE FIRSTMODE EXPERIENCES 
ANOMALOUS GVD, WHILE THE SECOND MODE PROPAGATES 
IN THE NORMAL DBPERSION REGIME. THE FIBER 
PARAMETERS ARE -0.06ps 2 /m,/? 1 ± 1 =0.06ps 2 /m. OTHER 
PARAMETERS ARE IDENTICAL TO THOSE LB ED FOR FIG. 1. 




With the increase of the parameter <5, it can be realized 
that the number of the spectral peak decreases and we 
observe only two sidebands. For the large values of 5, 
the bandwidth of the second sideband increases, while 
the peak gain of the first sideband decreases (FIG. 
3(a2)- (a4)). 




FIG. 3(a2) THE MI GAIN SPECTRA WITHOUT THE FWM 
NONLINEARITIES . THE FIRSTMODE PROPAGATES IN THE 
NORMAL DBPERSION REGIME, WHILE THESECONDMODE 
EXPERIENCES ANOMALOUS GVD. THE FIBERPARAMETERS 
ARE = 0Mps 2 /m,p± = - 0.06ps 2 /m. OTHER PARAMETERS 
ARE IDENTICAL TO THOSE USED FOR FE. 1. 




FIG. 3(a3) THE MI GAIN SPECTRA WITHOUT THE FWM 
NONLINEARITIES. THE FIRST MODE PROPAGATES IN THE 
NORMAL DBPERSDN REGIME, WHILE THE SECOND MODE 
EXPERIENCES ANOMALOUS GVD. THE FIBER PARAMETERS 
ARE 0.06ps 2 /m, /?* = - 0.06ps 2 /m. OTHER PARAMETERS 
ARE DENTIC AL TO THOS EUSED FOR FE . 1 . 




FIG. 3(al) THEMI GATN SPECTRA WITHOUT THE FWM 
NONLrNEARITiES. THE FIRSTMODE PROPAGATES IN THE 
NORMAL DBPERSION REGIME, WHILE THESECONDMODE 
EXPERIENCES ANOMALOUS GVD. THE FIBERPARAMETERS 
ARE = 0.06ps 2 /m, /?*= - 0.06ps 2 /m. OTHER PARAMETERS 
ARE IDENTICAL TO THOSE USED FOR FE. 1 



FIG. 3(a4) THE MI GAIN SPECTRA WITHOUT THE FWM 
NONLTNEARrnES. THE FIRST MODE PROPAGATES IN THE 
NORMAL DBPERSDN REGIME, WHILE THE SECOND MODE 
EXPERIENCES ANOMALOUS GVD. THE FIBER PARAMETERS 
ARE 0.06ps 2 /m, /?£=- 0.06ps 2 /m. OTHER PARAMETERS 
ARE DENTIC AL TO THOSE USED FOR FE. 1. 
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Furthermore, Fig. 3 (bl)-Fig. 3 (b4) show also that with 
the presence of the FWM coefficients, the spectral 
peaks are different with those obtained in Fig. 3 (a). 




FIG. 3(bl) THE MI GAIN SPECTRA IN THE CASE WHERE THE FWM 
NONLrNEARITIES ARE TAKEN INTO ACCOUNT. THE FIRST 
MODE PROPAGATES IN THE NORMAL DISPERSDN REGIME, 
WHILE THE SECOND MODE EXPERIENCES A NOMA LO LB GVD. 
THE FIBER PARAMETERS ARE /?* = 0.06ps 2 /m, /?±= - 0.06ps 2 /m. 
OTHER PARAMETERS ARE IDENTICAL TO THOSE LB ED FOR 
FIG. 1 




FIG. 3(b2) THE MI GAM SPECTRA IN THE CASE WHERE THE FWM 
NONLrNEARITIES ARE TAKEN INTO ACCOUNT. THE FIRST 
MODE PROPAGATES IN THE NORMAL DISPERSDN REGME, 
WHILE THE SECOND MODE EXPERIENCES A NOMA LO LB GVD. 
THE FIBER PARAMETERS ARE /?* = 0.06ps 2 /m, /?*= - 0.06ps 2 /m. 
OTHER PARAMETERS ARE IDENTICAL TO THOSE LB ED FOR 
FIG. 1 




FIG.3(b3)THE MI GAIN SPECTRA IN THE CASE WHERE THE FWM 
NONLrNEARITIES ARE TAKEN INTO ACCOUNT. THE FIRST 
MODE PROPAGATES IN THE NORMAL DISPERSDN REGME, 
WHILE THE SECOND MODE EXPERIENCES A NOMA LO LB GVD. 
THE FIBER PARAMETERS ARE = 0.06ps 2 /m, = - 0.06ps 2 /m. 
OTHER PARAMETERS ARE IDENTICAL TO THOSE LB ED FOR 
FIG. 1 




FIG. 3(b4) THE MIGAIN SPECTRA IN THE CASE WHERE THE 
FWM NONLINEARITiES ARE TAKEN INTO ACCOUNT. THE 
FIRST MODE PROPAGATES IN THE NORMAL DBPERSION 
REGIME, WHILE THE SECOND MODE EXPERIENCES 
ANOMALOUS GVD. THE FIBER PARAMETERS ARE /?* = 
0.06p> 2 /m, /?±= - 0.06ps 2 /m. OTHER PARAMETERS ARE 
IDENTICAL TO THOSE LB ED FOR FIG. 1 

Numerical Simulations in Different GVD 
Domains 

Direct numerical simulations were run in order to 
understand the dynamics of CW states under the MI 
in the nonlinear regime. We have solved Eq. (1) using 
a split-step Fourier method applying the fast Fourier 
transforms [24]. The numerical simulations have been 
carried out using various step sizes down to 0.0001 
throughout the z direction, with up to 1024 points 
along the t direction. The CW state was initially 
perturbed as: 



x oi = y^"(l±0.005)cos(nt) 
Xfi = V^(l± 0.005)cos(flt) 



(8) 



whereP = 1.5 is the power of the CW and Q = 0.2 is 
the angular frequency of the weak sinusoidal 
modulation imposed on the CW. Examples of different 
numerical results are shown below. 

Casel: ^ < O./fo < 0,/^ > O.jSfi > 

FIG. (4) shows the evolution of the CW solution given 
by Eq. (8). In FIG. 4 (a), we have neglected the FWM 
nonlinearities and we can see here a formation of a 
periodic wave train. The period of these pulses is the 
same as that of the initial perturbation and equal to 
2ti/Q. 
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FIG. 4 (a) DIFFERENT PROPAGATION RESULTS WITHOUT THE 
FWM NONLINEARITIES . THE FIBER PARAMETERS ARE /?± = 
-0.06ps 2 /m, /?± = 0.06ps 2 /m. OTHER PARAMETERS ARE 
IDENTICAL TO THOSE USED FOR FIG. 1 

A similar result had been obtained in a system of two 
parallel- coupled fibers modeled by two linearly 
coupled complex Ginzburg-Landau equations [25] and 
in a system of higher -order complex Ginzburg-Landau 
equation, with the fourth-order dispersion and cubic- 
quintic nonlinear terms [26]. But after a certain 
distance of propagation, the wave trains become 
unstable. By comparison, in FIG. 4(b) where the FWM 
coefficients are taken into account, we observe that the 
pulse can still propagate stably without any distortion. 
But at the end of the propagation, the number of the 
pulses is increased by a factor two. 




FIG. 4(b) DIFFERENT PROPAGATION RESULTS WITH THE FWM 
NONLINEARITIES. THE FIBER PARAMETERS ARE jS* = 
-0.06ps 2 /nvft ± 1 = 0.06p> 2 /m. OTHER PARAMETERS ARE 
IDENTICAL TO THOSE USED FOR FIG. 1. 



Case 2: ftj > 0,0^ > 0,/^ < 0,/?^ < 

In FIG. 5 (a) where we have neglected the FWM 
nonlinearities, we show an example of the first 
situation where a periodic wave trains is formed. But 
as in the previous case the pulses are distorted at the 
end of the propagation. 




FIG. 5 (a) DIFFERENT PROPAGATION RESULTS WITHOUT THE 
FWM NONLINEARITIES. THE FIBER PARAMETERS ARE/?± = 
0.06ps 2 /m, /?*= - 0.06ps 2 /m. OTHER PARAMETERS ARE 
IDENTICAL TO THOSE USED FOR FIG. 1 



Another example of the evolution is shown in FIG. 5(b) 
where the FWM nonlinearities are taken into account. 
We can see the amplitude of the wave trains increases 
progressively with the propagation distance. However, 
contrary to the previous case, the number of the pulses 
is constant. 




FIG. 5(b) DIFFERENT PROPAGATION RESULTS WITH THE FWM 
NONLINEARITIES. THE FIBER PARAMETERS ARE/?± = 0.06ps 2 /m, 
/?*= - 0.06ps 2 /m. OTHER PARAMETERS ARE IDENTICAL TO 
THOSE LB ED FOR FIG. 1 



1.24 




t 



FIG. 6(a) DIFFERENT PROPAGATION RESULTS WITHOUT THE 
FWM NONLTNEARTHES. THE FBERPARAMETERS ARE/?* = 
-0.06ps 2 /m, /?*= - 0.06ps 2 /m. OTHER PARAMETERS ARE 
IDENTICAL TO THOSE LB ED FOR FIG. 1 
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Case 3: < 0,/? "i < 0,/?^ < 0,^ < 

The obtained results are completely different from 
those obtained in previous case. For example, when 
we neglect the FWM nonlinearities, we observe in FIG. 
6 (a) that at the beginning of the propagation, the 
waves train are stable. But, at the end of the 
propagation, these wave trains are distorted. 

By taking into account the effect of the FWM 
nonlinearities in FIG. 6 (b), we observe that the waves 
train still propagate stably without any distortion. 




1.2 5Q 



FIG. 6(a) DIFFERENT PROPAGATION RESULTS IN THE 
PRESENCE OF THE FWM NONLINEARITIES. THE FIBER 
PARAMETERS ARE /?± = -0.06ps 2 /m, 0*= - 0.06ps 2 /m. OTHER 
PARAMETERS ARE IDENTICAL TO THOSE USED FOR FIG. 1 

Conclusions 

The full treatment of MI of continuous waves for a 
system of coherently CNLSE relevant to the 
polarization of light beams has been studied. In 
contrast to the previous works, the finite birefringence 
effect and phase dependent nonlinearity (i.e. the FWM 
terms) have been considered. The results indicate that 
both the GVM and the FWM terms will alter the gain 
spectra of MI considerably in terms of shape, position 
and especially the number of the spectral regions, but 
in different ways. Moreover, the nonlinear evolution 
of MI for the presence of phase dependent 
nonlinearity has been studied numerically using split- 
step Fourier method. The most notably change occurs 
when both beams propagate in the anomalous 
dispersion regime. We have observed here that when 
the FWM nonlinearities are omitted, the MI leads to a 
pattern in the form of periodic pulses array which are 
distorted at the end of the propagation. However, 
when the FWM nonlinearities are added, we have 
obtained that periodic pulses might propagate without 
any distortion. In order cases in which two beams 
propagate in the anomalous dispersion regime and the 
two others in the normal dispersion regime, the 
situation is quite different. The undistorted pulses are 



observed only in the case where the FWM 
nonlinearities are omitted. 
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APPENDIX 

Matrix elements of Eq. (7) are given by 

m ll = ~ m 22 = -™44 = m 33 = -Pol ^ + ff 0i n_ 
O 1^01' m 12 = ™34 = ™23 = -™21 = -™43 = -01*01' 



H2 



02,1*01*01' 



- (°3 A 1*01*11 + 06.1*01*11 



m 17 = m 1H = m, q = m qfi = - m 77 = - m 7R = -m AK = 



-m. 



L 35 



L 36 



27 



'28 



L 45 



(0 4 >oT*iT + 0JP*P* ), 



l 84 



(03j*S*£ + Q 6 Jp£p?i )> 



m 53 = m^ A = -m fi , = - m fii = m 71 = m 77 = -m H1 = 



L 54 



L 63 



64 



'71 



l 72 



~m 82 = -(04/11 *£ + 05j* , ii* 5 ol )i 

^55 = -^66 = ^77 = ^88 = ~Pu^ 2 + ~ 01*11 ' 

— - i 78 = — m 87 = —Q\P{{, 



l 56 



l 65 



l 57 



l 67 



l 68 



l 75 



l 76 



l 85 



'"■b - 02.1*11*11' 
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